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Abstract: We show explieitly, by using astrophysieal data plus reasonable assumptions 
for the bulk viscosity in the cosmic fluid, how the magnitude of this viseosity may be 
high enough to drive the fluid from its position in the quintessence region at present time 
/ = 0 aeross the barrier w = —1 into the phantom region in the late universe. The phantom 
barrier is aceordingly not a sharp mathematieal divide, but rather a fuzzy eoneept. We also 
ealeulate the limiting forms of various thermodynamical quantities, including the rate of 
entropy produetion, for a dark energy fluid near the future Big Rip singularity. 
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1. Introduction 

Recent years have seen an increased interest in viscous cosmology. This contrasts the traditional 
approach, in which the cosmic fluid has usually be taken to be ideal (nonviscous). Erom a 
hydrodynamicist’s point of view this is to some extent surprising, since there are several situations 
in fluid mechanics, even in homogeneous space without boundaries, where the inclusion of viscosity 
becomes mandatory. Especially, this is so in connection with turbulence phenomena. As is known, to 
the first order in deviation from thermal equilibrium there are two such coefficients, the shear viscosity 
r], most often being the dominant one, and the bulk viscosity C [!]• Because of the assumption about 
spatial isotropy in cosmology, the shear coefficient is usually omitted. (As a parenthetical remark we 
note that this first-order theory away from thermal equilibrium means in effect acceptance of the Eckart 
1940 theory [2], implying noncausality. By taking into account second order deviations from thermal 
equilibrium one can, however, obtain a causal theory respecting special relativity. Pioneering articles on 
causal fluid mechanics are those of Muller [3], Israel [4], and Israel and Stewart [5].) 
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Viscosity coefficients have so far been considered before, both for the past {t < 0) and for the 
future {t > 0) universe. Sinee a large amount of knowledge is available about the past universe, it 
is natural to apply methods from partiele physios in order to understand the visoosity’s influence upon 
the evolution. Thus, Hogeveen et al. [6] give results for both shear and bulk viseosity in the very 
early universe, oaloulated from kinetio theory. The general tendeney seems to be that the modifieations 
from viseosity to the development of the early universe are moderate or small. In the future universe, 
however, as we will focus on below, muoh less is known, and it beoomes therefore natural to make use 
of phenomenological suggestions for the magnitude of the viseosity. Some artieles and reviews dealing 
with viseous eosmology from a phenomenologieal point of view are Refs. [7-12]. 

We will in the following eonsider the future evolution of the eosmie fluid, where for simplieity this is 
taken to eonsist of one eomponent only. This is the dark energy fluid endowed with an equation-of-state 
parameter w. Again for simplicity, and also because of the desirability to keep the number of input 
parameters in the theory restrieted, we shall assume that te is a eonstant. Thus, taking the equation of 
state to be homogeneous, we ean write 


p = wp, w = —1 + a. (1) 

Astronomical observations show that w is close to unity, thus, a very small. Aeeording to the 2015 
Planek data, w = —TOlO^o ggo [13] (see Table 5). The region —1 < te < —1/3 is the quintessence 
region, while te < — 1 is the phantom region. The ease w = —1 is ealled the phantom barrier. 

In this paper we will deal with the phantom region and the phantom barrier. For some years it has 
been known that if the eosmie fluid starts out from some value of w lying in the phantom region, it 
will eneounter some form of singularity in the remote future. The most dramatie event is ealled the Big 
Rip, in which the fluid enters into a singularity after a finite time span [14-16]. There are also softer 
variants of the future singularity where the singularity is not reaehed until an infinite time, ealled the 
Little Rip [17-19], the Pseudo Rip [20], and the Quasi Rip [21]. 

In view of the mentioned smallness of |a| it becomes natural to ask: if the eosmie fluid starts from 
present time f = 0 in the quintessenee region {a > 0), has it any possibility to slide through the phantom 
barrier into the phantom era and thereafter inevitably be drawn into the future singularity? This problem 
was analyzed already in 2005 [22] with the following result: if the magnitude of the bulk viseosity is large 
enough, and if the viseosity varies with the energy in a way that is natural physieally, sueh a transition 
is aetually possible. The analysis in Ref. [22] was however restrieted to qualitative eonsiderations, as 
quantitative information about the viseosity was not available. Reeently this situation has improved, due 
to the analysis of Wang and Meng [23] and others on the influenee from viscosity on the Friedmannn 
equations giving the evolution of the Hubble parameter, H = H{z), z being the redshift. Introduction 
of the quantitative element is the main point in the present study. As will be shown, reasonable values 
of the viseosity inferred from the experiments are just of the right order of magnitude neeessary to make 
the transition through the phantom barrier possible. 

It has reeently eome to our attention that there exists a reeent series of papers by Velten et al. dealing 
with viseous cosmology [24-27], of definite interest for our work. We shall have the opportunity to 
compare with some of those results later. See also the related paper of Nunes and Pavon [28]. 
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In the next seetion we present the general formalism, list the various aetual forms for the relationship 
C = C(p)> diseuss the fluid’s possibilities for phantom barrier passage. In Seetion 3 we eonsider 
further eomparisons with experiments. In Seetion 4 we diseuss, under simplifying assumptions, limiting 
forms for various thermodynamie quantities near the future singularity. 


2. Basic Formalism. Various Assumptions for the Bulk Viscosity 

We first have to sketeh the basie formalism as outlined, for instanee, in Refs. [22] and [29]. We start 
from the standard FRW metrie 

ds^ = —dt^ + a^{t)dy^, ( 2 ) 

when the spatial eurvature k is set equal to zero. 

The Friedmann equations, with 6 = 3H = 3d/a the sealar expansion, are 


= 247rG' p + 


A 


8tiG 


0 + 1:0^ = -UttG { p - C(p)0- 


A 

’ 


and the energy eonservation equation is 


p + {p + p)0 = C{p)0‘^- 


(3) 

(4) 

(5) 


We set heneeforth the eosmologieal eonstant A = 0. The governing equation for the sealar expansion 
then beeomes 

0 + ^a02 - 127rG'C(p)0 = 0, (6) 

whieh ean in view of Equation (3) be rewritten as 

p + \/247rGap^/2 _ 2A'kGC,{p)p = 0. (7) 

The solution is 


t = 


1 


dp 


v247rG J pQ p3/2 


-a + \/24:7rGC{p)/y/p 


( 8 ) 


where we integrate from present time t = 0 (with density po), into the future. Note that the sign of a as 
defined in Equation (1) is the opposite of that used in Ref. [22]. 

We will now analyze different options for C(p)- 


2.1. ( Equal to a Constant 

In this ease it is most eonvenient to go baek to Equation (6) and solve direetly for 0: 

006*/*'^ 


0(f) = 


1 + \a0Qtc.{e^Ac - 1 )’ 


(9) 


where tc = {12'kGC,) ^ is the “viseosity time”. 
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Correspondingly, we obtain for the seale faetor 


Q,(t) — Q,q 


1 + -a6otc(e 


t/tc 


1 ) 


2 

3a 


( 10 ) 


and for the density 


pit) = 


Poe 


2t/t^ 


[l + - 1)]' 


( 11 ) 


We thus see that if the universe starts out from the phantom region (a < 0), it will inevitably be driven 
into the Big Rip singularity at a well-defined rip time t = tg, where 


tg = tjn(l+ ^ ) . (12) 

Both 9{t),a{t), and p{t) diverge at this instant. However, if it starts out in the quintessenee region 
(a > 0), the behavior is quite different: we see that in the remote future 

2 

0(f) —y —, t —y oo, (13) 

ate 

a{t) —)■ oo, f —)■ oo, (14) 

pit) t^oo. (15) 

ia9otc)^ 

Observe that the finite limiting values for 9{t) and p{t) are viseosity dependent. For an ideal fluid 
{tc —)■ oo), these two quantities go to zero when t ^ oo. 

We see that in this model the ease a = 0 aets as a sharp boundary; the behavior of the fluid when a is 
positive is very different from the behavior when a is negative. This indieates that the assumption about 
eonstant ( is over-idealized; a real physieal fluid is not expeeted to behave in sueh a singular way, at least 
not so in a elassieal theory. We therefore turn to models whieh are physieally more plausible. There are 
two models of that kind, to be eonsidered sueeessively below. 


2.2. ( Proportional to 9 

In the late universe, espeeially near the future singularity where the fluid motion is vigorous and 
the oeeurrenee of turbulenee highly likely, it is natural to assume that the bulk viseosity inereases with 
respeet to time. We eonsider first a model where C oc 

Cip) = ti9 = ri^247rGp, (16) 


with Ti a eonstant. Then Equation (8) yields 

1 1 
y/f^ y/P 

with pq referring to the present time t = 0. If the universe starts out from the phantom region, the 
development into a Big Rip singularity is inevitable. But if it starts out from the quintessenee region, the 



-J_^_( 

^/2A^^G -Oi + 2 AtiGti \ 
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development depends critieally on how large the coefficient ti is. The condition for passing through the 
phantom barrier and into the phantom region is seen to be 


247rGri 


> a, 


(18) 


when expressed in dimensional units. If this condition holds, the Big Rip (p = oo) is reached in a 
finite time. 

This possibility, pointed out already in Ref. [22], can now be supplemented by data inferred from 
experiments, though only on the order-of-magnitude level. We shall first consider the analysis of Wang 
and Meng [23]. They include various assumptions for the bulk viscosity in the early universe in the 
Friedmann equations and compare the theoretical curve for H = H{z) with a number of observations. 
The comparison gets somewhat complicated, due to the various parameters involved, and also because 
of allowance of a more general equation of state than the form w = — 1 + a which we have chosen above 
for simplicity. We intend to analyze this issue more closely in a forthcoming paper [30] . For our present 
purpose the relevant information can however be drawn relatively easily: a central quantity is B, which 
we shall define as 

127rG'Co 

B = - (19) 

Here UttG/c^ = 2.79 X 10 m/kg, and the dimension of i? is [i?] = s Again, in the subscript 
zero refers to present time / = 0. 

From Refs. [23] and [30] we infer that, given our simple form for the equation of state, the best 
agreement with observational data is obtained when 


B 


rsj 


50 


km 

sMpc 


1.62 X 10"^V^ 


( 20 ) 


[Note: there is an apparent difficulty here because the values of B turn out to be negative. This may be 
due to two causes; either that (i) the bulk viscosity introduced in Friedmann’s equations is of a formal, 
rather than a physical, nature; (ii) or that this behavior is related to our choice for equation of state. We 
think the latter option is true here. The analysis turns out to be quite sensitive to the different equations 
of state one adopts for the fluid. More complicated forms for the equation of state easily lead to positive 
values for B. In any case, the magnitudes for |i/| turn out to be roughly in agreement with Equation (20). 
This is sufficient for our purpose.] 

From Equations (19) and (20) we estimate, on basis of the given information. 


Co ~ 5 X 10^ Pa s. 


( 21 ) 


This value is probably somewhat high. We may compare with the formula given by Weinberg [7] for the 
bulk viscosity in a photonic fluid. 


C = daradP'^T/ 



( 22 ) 


where Orad = 7r^k^/15h^c^ is the radiation constant and Tf the mean free time. Considering the 
microwave background at 3 K slightly out of equilibrium with pressure-free matter and estimating the 
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free time to be the inverse Hubble radius, Tf = we obtain ( ~ 10^ Pa s. Let us in the following 

simply assume, as a eonservative estimate, that Co lies in the region 


10^ Pa s < Co < 10^ Pa s. 

(23) 

Now go back to Equation (18) and insert for ti present-time quantities, Ti 
condition for phantom barrier passage can then be written 

= Co/^o> with 9q = SHq. The 


(24) 

With Hq = 67.80 km s“^ Mpc“^ = 2.20 x 10“^® we obtain 


Co > (1.18 X 10^) a. 

(25) 


We gave above (after Equation (1) the value of w aeeording to the 2015 Planek data; this yields the 
maximum value to be tUmax = —1.019 + 0.075 = —0.944 or Omax = 0.056. In this ease, the eondition 
for driving the fluid aeross the phantom barrier is 

z47rG 

This leads us to the important eonelusion that the the estimates for the present viseosity given in 
Equation (23) are roughly high enough to drive the fluid through the phantom barrier, even in the most 
extreme quintessential ease w = W max . 

2.3. C Proportional to 9^ 

Another option would be to assume that the bulk viseosity varies with 6 as 

C(p) = T2 6"^ = T2- 2AttGp, (27) 

with T 2 a new eonstant. This amounts physieally to attaehing more weight to the stages where the fluid 
density is high. In this ease we obtain, instead of Equation (17) [31], 

-2 dx 

\/2AtiG J./p^ [a - (247rG)3/2r2 x] 

2 ,1/1 _ 1_\ mnG)^ f^o a-(247rGr/^T2^ \ 

\/24^V VVP a-{247rG)^/^T2^J^' 

Again, we find that the future development of the universe is eritieally dependent on the value of the 
equation-of-state parameter, assumed eonstant in the present model. If the initial stage (t = 0) admits 
one to put a < 0, then a future Big Rip singularity will be eneountered after a finite time 

~ H V24vrG>o' 

However, if the universe starts from the quintessenee region at f = 0, the expression (27) shows a 
eomplieated behavior indieating that the expression as sueh loses physieal interest. The reason for this 
is that the assumption w = —1 + a that we have made in this paper for simplieity is incompatible with 
the viscosity ansatz (27) above. One may easily change the equation of state into a more general form 
w = w{p), and in that way maintain the same sliding property of the cosmic fluid [31]. 
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3. Further Comparisons 


In our numerical considerations above, we focused on the analysis of Wang and Meng [23]. It is 
desirable to compare with other investigations also. As mentioned at the end of Section 1 there are 
several recent papers in this direction. A paper of great interest is that of Velten et al. on dark matter 
dissipation [24]. A main outcome of their analysis of current data from supernovae, baryon acoustic 
oscillations and cosmic microwave background, was that dark matter is allowed to have a bulk viscosity 
< 10^ Pa s. This was found to be consistent also with the integrated Sachs-Wolfe effect (which 
sometimes is a difficult point for viscous cosmology theories). The agreement with our estimate in 
Equation (23) is actually seen to be quite good. 

It should be mentioned here for the sake of comparison that the nondimensional quantity ^ in Ref. [24] 
is defined to be essentially the same as our B/Hq, 


i 


2B 


The authors represent the bulk viscosity in the form 


c 



(30) 


(31) 


with emphasis on the cases o = 0 and z/ = —1/2. 

Another noteworthy point in connection with Ref. [24] is that the viscous theory was found to be able 
to form galactic dark halos (within the hierarchial scenario) only if ^ <^0.2. This means 


Co < 7.9 X 10® Pa s. 


(32) 


This behavior should accordingly be quite compatible with the range given in Equation (23). Most of the 
cases shown in Eigure 2 in Ref. [24] are however given for lower viscosities. As an example, the typical 
value of C = 2 X 10“'^ given for galactic clusters corresponds to Co = 7.9 x 10^ Pa s. A general tendency 
is that the higher viscosities belong to the larger scales (galaxy clusters). 

Similar results are obtained in Ref. [27] on structure formation in viscous universes. Eor instance, if 
dark matter has a viscosity of C = 10“®, a strong growth suppression is observed at small scales. 


4. Remark on Entropy Production. Behavior of Thermodynamic Quantities near the Future 
Singularity. 


We may start by applying the standard formula for the rate of entropy produced per particle (subscript 
zero again referring to present time). 


do 


Co q2 
noksTo ° 


(33) 


We may interpret the universe as a “gas”, letting a “particle” be a galaxy, and we may insert for the 
number density of galaxies the value [32] 


n ~ 10 m 


(34) 
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Inserting a typical bulk viscosity of Co ~ 10^ Pa s, and using 6*o = 3i/o = 6.60 x 10“^® together with 
To = 3 K, we then estimate from Equation (33) an enormous value for the entropy growth which escapes 
a natural physical interpretation. This may be related to the difficulty of assigning a viscosity to a gas 
of galaxies. We find it better to leave this point, and analyze instead the behavior of the thermodynamic 
quantities near the future Big Rip, which is a topiof definite physical interest. Let us assume that the 
universe starts out with a value of a that is small, but nevertheless negative. For simplicity we will limit 
ourselves to the case of constant viscosity, as considered in Section 2.1. Thus, we assume henceforth 

w = —1 + a = constant, C = constant. (35) 


From Equations (9)-(l 1) we obtain 


6{t) -)■ 


2 1 
Ci\ts — t' 


t ^ tg, 


ci{t) —y 


ta 


1 + \ \Ot\9Qtc 
4 


2 

3|c«l 


Oq 


(t, -t)2/3|«l’ 


t —)■ tg, 


p{t) 


Po 


{aeoY {tg-ty^ 


t —>■ tg] 


(recall Equation (12) for the rip time tg). 

We need also the corresponding limit for n{t). From the continuity equation {nU^).^ 
the local rest frame fi = —n6, which implies 


n{t) —)■ no 


+ ^\a\9otc^ 1^ 


t)l^, t^tg. 


(36) 

(37) 

(38) 


0, we have in 


(39) 


This expression actually goes to zero, at the singularity. 

Finally, we shall find the time dependence of T{t) near the singularity, under the assumption that the 
viscosity is “small”. This is taken to mean that the approximation 


a\9otc 2> 1 


(40) 


holds. We may actually check this in the present case, by taking ( = 10^ Pa s as a typical value. In 
dimensional units we find 

Thus, as long as |a| is of order 10“^ or less, the approximation (40) holds reasonably well. In view of 
Equation (12), this means that 

tg <^tc. (42) 


We can then find T(t) from the standard formula for nonviscous flow, using |q!| <C 1, 


T{t) = To - 
ao 


3(l«l+l) 


To — 


Tn 


(X 


t ^ tg. 


aoj (t.-f)2/H’ 

Making use of these expressions, we find for the rate of entropy production per galaxy 


a{t) = 


C 


n{t)kBT{t) 


9\t) 


C 


(43) 


«2 (tg-tr 


t ^ tg. 


(44) 
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The entropy produetion thus diverges at the Big Rip, as we might expeet. A noteworthy property of this 
expression is that |a|, though present in the prefaetor, is absent from the exponent. We see that a{t) 
behaves similarly as the density p{t). Equation (38), near the singularity. 

5. Summary 

A main lesson from this investigation is that the phantom barrier w = —1 should not be regarded as a 
sharp mathematieal divide. The eosmie fluid is a eomplieated physieal system, most probably endowed 
with properties sueh as viseosities sueh as eommon elsewhere in fluid meehanies. Beeause of spatial 
symmetry, we have ignored the shear viseosity and maintained the bulk viseosity only. As we have seen 
in Seetion 2, the estimate (23) for the present-time bulk viseosity makes a passage through the phantom 
barrier quite possible in the late universe, even if it starts at present time t = 0 from the quintessenee 
region te > —1 (« > 0), at a maximum value w = tUmax (= 0.056) inferred from the Planek data 
(2015) [13]. Our eomparison with experimental results was based mainly on the analysis of Wang and 
Meng [23], although eomparison with other analyses of Velten et al. [24,27] gave eomparable results. 
This agreement is eneouraging, and indieates that our numerieal estimates are on the right traek. 

The qualitative derivations of the limiting forms of thermodynamie quantities near the future 
singularity t = tg in Seetion 4 were made on basis of the simplifying assumptions that both w and 
the bulk viseosity C are eonstants. Near the singularity Equation (44) shows that the entropy produetion 
diverges, as expeeted, whereas it eould hardly have been seen beforehand that the exponent turns out to 
be independent of |q!|. 

It should be emphasized that our theory is restrieted to a one.eomponent model. This is motivated, 
of eourse, by the dominant role of the dark fluid eomponent. A more aceurate deseription would be 
obtained by assuming n eomponents of the eosmie fluid, 

n n 

P = ^Ph P = ^ ^iPi^ (45) 

i=l i=l 

{P being the total pressure, Wi =eonstants), we ean insert these into the eontinuity equation 

adap{a) + 3[p{a) + P] = 3(6, (46) 

here expressed in terms of a. We intend to study this model more elosely in Ref. [30]. 

As a elosing remark, it eould be mentioned that the topie we have eonsidered in this paper is related 
to the so-ealled quintom eosmology. The simplest quintom model is aehieved by introdueing two sealar 
fields, one of them being quintessenee and the other a phantom field. This kind of formulation arises 
from the desire to find alternative, and perhaps better deseriptions, of the universe than inflationary 
theory. Crossing of the phantom barrier w = —1 is an important ingredient of the theory, and ean 
even be deseribed within the framework of braneworld eosmology. An extensive review of quintom 
eosmology is given by Cai et al. [33]. A reeent paper relating quintom matter to the so-ealled emergent 
gravity ean be found in Ref. [34] . 
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